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ABSTRACT 

We report turbulence effects on magnetic reconnection in relativistic plasmas using 3-dimensional relativistic 
resistive magnetohydrodynamics simulations. We found reconnection rate became independent of the plasma 
resistivity due to turbulence effects similarly to non-relativistic cases. We also found compressible turbulence 
effects modified the turbulent reconnection rate predicted in non-relativistic incompressible plasmas; The re¬ 
connection rate saturates and even decays as the injected velocity approaches to the Alfven velocity. Our 
results indicate the compressibility cannot be neglected when compressible component becomes about half of 
incompressible mode occurring when the Alfven Mach number reaches about 0.3. The obtained maximum 
reconnection rate is around 0.05 to 0.1, which will be able to reach around 0.1 to 0.2 if injection scales are 
comparable to the sheet length. 

Subject headings: magnetic fields, magnetohydrodynamics (MHD), magnetic reconnection, relativistic pro¬ 
cesses, plasmas, turbulence 


1. INTRODUCTION 


Magnetic reconnection is known as a process responsible 
for a very efficient magnetic field dissipation in many plasma 
phenomena. In particular, it is expected to play an important 
role for the acceleration of relativistic outflow in high energy 


plasmas, such 

as relativistic jets (Blandford & Znajek[[1977 

Begelman et al. 1984 Komissarov et al. |2007t Kino et al. 

2015|l, pulsar \ 

& Kirk||2001t 

vind (Kennel & Coroniti 1984a|b| Lyubarsky 

Kirk & Skjseraasenj 20()3|l, and gamma-ray 

bursts (GRB) 

Lyutikov & Blandford[[2003[ [Zhang & Yan 


tion ( |Sweet|1958t|Parker|1957[) predicts that magnetic recon- 
nection becomes very slow in high magnetic Reynolds num¬ 
ber plasmas (/?„! ~ 10'®), and fails to explain observed dissi¬ 
pation timescale in space and astrophysical phenomena. To 
solve this problem, a lot of efforts have gone into finding a 
fast-reconnection process that does not depend on the value 
of resistivity. Using the equation of continuity, the reconnec¬ 
tion rate can be expressed as 


Oa Pin Oa L 


( 1 ) 


where the subscript “in” and “s” indicate the inflow and out¬ 
flow region, respectively, Vin, Vs are the inflow and outflow ve¬ 
locity, respectively, ca is the Alfven velocity, p is the mass 
density, 6 is the sheet thickness, and L is the sheet length. 
This equation shows that fast reconnection processes can be 
obtained by increasing the density ratio: Ps/Pin, (Brunei 
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tio of 

sheets: d/L ( 

Biskamp|1986l [Shibata & Tanuma|2001 

Loureiro et al.|2007 

[Bhattacharjee et al.|2UU91 Uzdensky et al.j 

2010l[Takamoto|20l3 

Sironi & Spitkovsky|2014 1 . 


accelerate magnetic held annihilatio n ( Matthaeus & Lamkin 
1985 |Eyink|[20TT] [Takamoto et all2012[ [Higashimori et al. 


2013 I. In particular, many astrophysical objects are consid- 
ered to be high Reynolds number plasma, and it is natural to 
assume those plasma are in a turbulent state Q It was theoreti¬ 
cally suggested that strong Alfvenic turbulence also increases 
the sheet aspect ratio, and the reconnection rate becom es inde¬ 
pendent of the resistivity ( [Lazarian & Vishniac| ( |1999] l, hence¬ 
forth LV99). LV99 predicts the following expression of recon¬ 
nection rate: 



( 2 ) 


where I and v; are the energy injection scale and velocity dis¬ 
persion of turbulence at the injection scale, respectively. This 
was ex amined using magn etohydrodynamics (MHD) simu¬ 
lation ( [Kowal et al.[[200^ . However, the numerical work 
was limited only in the non-relativistic incompressible regime 
with plasma (3 larger than unity, and its applicability to rel¬ 
ativistic Poynting dominated plasma with relativistic turbu- 


* In Poynting dominated plasmas, relatively strong turbulence may be 
abl e to be induced by yarious ways, e.g., Kruskal-Schwarzchild type instabil- 
ity |Lyubarsky|2010| and Richtmyer-Meshkoy type instability |Inoue|2012| 
at a shock tront. wh ich induces turbulence with yelocity dispersion up to 
< 1.5/v^ (Takamoto et al.|2012|; th e tearing instability with tur¬ 
bulent yelocity Lorentz factor, 7ti,rb — \/q:(t/ 2, where a is the magnetization 
parameter defined later, a is the sound yelocity, and a is energy conversion 
factor from magnetic field into kinetic energy. 
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lence was unclear, which is very impo rtant in the context of 
high energy astrophysical phenomena (|Lyutikov & Lazarian 
|2fei[Kadowaki et al.|2015] i. 

In this paper, we extend the previous work to relativis¬ 
tic plasma including both matter and Poynting dominated 
plasma. We also investigate effects from compressibility on 
reconnection rate. In Section 2 we introduce the numerical 
setup and the method for the turbulence injection. The nu¬ 
merical result is presented in Section 3, and its theoretical ex¬ 
planation is presented in Section 4. Their implications for 
some high energy astrophysical phenomena are discussed in 
Section 5. Section 6 summarizes our conclusions. 

2. SIMULATION SETUP 

We modelled the evolution of a current sheet in a turbulent 
flow using 3-dimensional resistive relativistic magnetohydro¬ 
dynamics (RRMHD). The in itial current sheet is modelled by 

the relativistic Harris sheet dHohll 19661 IKirk & Skiaeraasen 

- --- 


as 


|2003| l whose magnetic field is expresse 

B = Botanh[z/A]ei+Bce3,, 


(3) 


where A is the half-thickness of the initial sheet, and Bq and 
Be are the reconnecting magnetic field and guide field compo¬ 
nent, respectively. The pressure inside of the sheet is assumed 
to satisfy the pressure balance, and the upstream pressure is 
determined by the magnetization parameter a = B^/Airphc^ 
where h=l+(T/ (T-!))(/? / pc^) is the specific enthalpy of rel¬ 
ativistic ideal gas with T = 4/3, and p,p,caie the gas pressure, 
mass density, and the light velocity, respectively. The initial 
temperature is assumed uniform, Q = k^Tjme^ = 1 , where 
k^pn are the Boltzmann constant and particle rest mass, re¬ 
spectively. 

The evolution of the plasma is calcul ated using a 3- 
dimension al RRMHD scheme developed by Takamoto & In- 
|oue[ ( |20lT] ) which solves the full RRMHD equations in a con¬ 
servative fashion using the constrained transport algorithm. 
This allows us to treat the mass density, momentum, energy, 
and divergence of magnetic field to be conserved within ma¬ 
chine round-off error. The resistivity, ry, was assumed to be 
constant, typically p /Lc = 10""*. We follo wed the similar sim¬ 
ulation setup used in |Kowal et al.| ( (2009| l. The numerical box 
is assumed [-L/2,L/2J x [0,LJ x [-L,L] where L = 20A ^ . 
Note that the z-direction size of the numerical box is twice 
larger than x,y-direction to reduce the influence by turbulence 
on the reconnection inflow around z-boundaries. We divided 
the numerical box into the homogeneous numerical cells with 
size; A = L/512. The timestep size is set as: Af = O.lA/c. 
We set the periodic boundary condition in y-direction and free 
boundary condition x and z-direction. 

In our mo del, we drive turb ulence using a similar method 
described by|Mac Low|(|1999]l. We add a divergence-free 3- 
velocity held, 6v, and an electric held determined consistently 
to the injected velocity[^at time intervals Afinj in a box region 
located around the current sheet: [- X [ 0 ,L] X [- 4 , 4 ] 
where 4,4 a scale length that is sufficiently larger than 
the injected turbulence eddy scale; Afinj is chosen to be 
shorter than the eddy turnover time at the injection scale / 
: Afinj = V^Tii'inj.o where Vinj,o = 0.15c is a typical injec¬ 
tion velocity in this study. We note that the dynamics of 


^ Note that A is the initial half-width of the sheet, and is a constant. 

^ First, we splitted the electric field as: E = Edissip — r X B. Then, the ideal 
part is replaced by SE = -ti,cw x B where fnew is obtained by the relativistic 
addition law of v and 5v. 



Figure 1. The profile of magnetic field lines in a turbulent sheet and the gas 
pressure (back plane) in the case of tr = 5. The color bar is the gas pressure 
in the unit of the magnetic pressure in the initial inflow region. 


the turbulence becomes insensitive to the injection time in¬ 
terval as long as the injection time interval is arou nd the 
eddy turnover time at the injection scale. Following (|Inoue| 
et al. [2011|| and ([Takamoto et al. 2012|), the velocity field 


is described as; p5v‘ = ^sin(k • i*-l-i/f) where 7 is the 
Lorentz factor of the injected velocity, i covers {x,y,z}, and 
(j)l is a random phase. The one-dimensional power spec¬ 
trum of the velocity field is assumed flat, k^P(k) cx k^. The 


initial-field perpendicular Fourier components kj_ = ykj + k^ 

are chosen in a shell extending from ki - Ak to k; + Ak where 
4L/27r = 16, AkL/27r = 3. Note that this scale size is a little 
larger than the initial sheet scale, and can be well-resolved 
by our present resolution. The parallel-field wave number 
k|| =kx is determined by k|| = k_LVinj/cA where Vi„j = 
is the root-mean-square velocity. Since the injected turbu¬ 
lent veloci ty does not immediately foll ow the critical balance 
condition ( jGoldreich & Sridhar||1995]l, th e turbulence at the 
injected scale is weak (jGaltier et al.|2()00|l, which transits into 
the strong turbulen ce around the sheet width scale because of 
the energy cascade (jPerez & Boldyrev|2008[ Verdini & Grap-j 


|pin|2012] [Meyrand et al.|2015| l. lire weak MHD turbulence 

cascades the wave energy only perpendicular to the magnetic 
field, the turbulence strength, x = ta/tnl — k±yx/k\\CA, grad¬ 
ually increases up to unity, which results in the strong turbu- 
lence in the sheet (j Verdini & Grappin||2012[ [Meyrand et al.j 
2015 I. Note that the injected velocity Vinj is different from 
the velocity at the injection scale v;, which was first pointed 
out in LV99. At the injection scale, the weak MHD tur¬ 
bulence theory gives the following the energy cascade rate; 
Cinj ^ v//rNL ~ v//|//^Ca where tnl ~ (/_l/vz)^/(/||/ca) is 
the distortion time of Alfven wave packets. Combining this 
to the injected power: v?jj/Afinj, we obtain 


^ /^hnj/l 2 2 


caI]_ 


(4) 


where the inequality resulted from the excitation of compres¬ 
sion modes. In other words, Vinj is related to the injected 
power by some external force or free energy; on the other 
hand, v/ is the velocity resulted from the energy cascade of 
the weak MHD turbulence. We observed this relation in our 
simulations, and assume this in the following (see also ([Lazar- [ 
[ian & Vishniac|1999[[Kowal et al.|2009[ l). 
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3. RESULTS 


Figure 1 is a snapshot of the gas pressure and magnetic field 
lines in the turbulent sheet in the case of highly magnetized 
case; cr = 5. Differently from the Sweet-Parker sheet, the sheet 
is highly stochastic due to the turbulence, which induces a lot 
of reconnection points in the sheet and drives a fast reconnec¬ 
tion process. Figures 2 are the observed reconnectio n rates 
vr which is measured using a method proposed by |Kowal| 
|et al.|(|2009|l (see Equation (13) in this paper); this allows us 
to measure the effective value of Ey/Bo in the 3-dimensional 
case, which provides us a reconnection inflow velocity less 
contaminated by turbulent flows than the direct measure of 
inflow velocity, . The top panel shows the reconnection 
rate with respect to the injected turbulent velocity in various 
kinds of magnetized plasmas. This shows that the turbulent 
reconnection rate shows 3 characteristic behaviors depending 
on the injected turbulent velocity Vinj/ ca: (1 ) increasing re¬ 
gion following LV99; (2) saturation region giving maximum 
rate; (3) decreasing region. When the injected turbulent ve¬ 
locity is sufficiently small, incompressible approximation can 
be applied, and the reconnection rate grows following Equa¬ 
tion Q. On the other hand, when injected velocity becomes 
comparable to the Alfven velocity, turbulence becomes com¬ 
pressible and the reconnection rate deviates from the incom¬ 
pressible theory. Interestingly, the injection velocity Vinj /ca 
at the maximum rate becomes smaller as the magnetization 
parameter increases. We will discuss the relation of this ten¬ 
dency to the compressible effects in the next section. Note 
that the error bar in the panel seems decreasing with a. We 
consider this is because the kinetic energy of turbulence be¬ 
comes smaller comparing with the magnetic field energy as 
the magnetization parameter a increases. The bottom panel 
of Eigure 2 is the reconnection rate with respect to the dif¬ 
ferent Lundquist number. It shows the reconnection rate is 
independent of the Lundquist number, and determined by the 
turbulent strength. Note that the obtained maximum recon¬ 
nection rate is very fast, vr/ca ~ 0.05, and eve n compara¬ 
ble to the relativistic Petschek reconnection rate ( |Lyubars^ 
|2005j l. This maximum reconnection rate also indicate that it 
will be possible to reach around 0.1 to 0.2 if injection scales 
are comparable to the sheet length as indicated by Equation 
( 2 ) 0 . 


Figure 3 is the guide field dependence of reconnection rate 
in the case of Avinj = 0.15 ca and ctr = 5 where ctr is the 
magnetization parameter determined by reconnecting mag¬ 
netic field component Bq. We fixed the reconnected mag¬ 
netic fiel d and added the gui de field component. As was re¬ 
ported by|Kowal et al.|(|2009)l, the reconnection rate becomes 
independent of the guide field strength even in a relativistic 
Poynting-dominated plasma. In the case of Be/B q = 1, the 
time scale necessary for reaching the steady state becomes 5 
times longer than no guide field case. 


4. THEORETICAL CONSIDERATIONS 
4.1. Sheet Density 

The obtained reconnection rate in Figure 2 shows an inter¬ 
esting behavior owing to compressibility which cannot be ex¬ 
plained by incompressible theory, Equation Q. In the follow- 

^ In the following, the reconnection rate obtained using the method by 
Kowal et al (2009) is identified by the reconnection inflow velocity Vjn, which 
becomes valid in a statistically steady state because |Ej| ~ Vi^Bo in this case. 

^ Unfortunately, the dependence on the injection scale / is hard to test 
because of the limited inertial range of turbulence in the simulation. 
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Figure 2. Observed reconnection rate in its steady state. Top: Reconnection 
rate with respect to the injected turbulent velocity. Bottom: Reconnection 
rate with respect to the Lundquist number: S = Lca/t). 


AVi,j = 0.15CA 



Bq/Bq 


Figure 3. Observed reconnection rate in its steady state with respect to the 
guide field. Note that the guide field is added to the fixed reconnection field, 
Bx, so that the total magnetization pai'ameter increases with the guide field. 


ing, we give an explanation for the saturation and depression 
of the reconnection rate in high Alfven Mach number regime. 
Equation ([T]i indicates the compressible effects can be divided 
into 2 parts: (1) the density ratio between sheet and inflow re¬ 
gion Ps/Pin, (2) decrease of the sheet width (5/L. Note that <5 is 
the actual sheet thickness determined by the turbulence which 
is different from the initial thickness A. 

We begin with discussing the density ratio. Figure 4 plots 
Ps /Pin with respect to the injected turbulence velocity in the 
matter and Poynting dominated cases cr = 0.04 and 5, respec¬ 
tively. They show that the density ratio decreases linearly with 
the turbulent strength. This can be understood from the con- 
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Figure 4. The density ratio between that of inflow and sheet region: ps/pin 
in the cases of Poynting dominated case: <7 = 5 and matter dominated case: 
cr = 0.04. The ratio decreases with increasing the turbulent strength due to 
the compressible effect. 


servation of energy flux; 

Pin^in^ (1 “I” Pin( 1 


(5) 


We assumed a non-relativistic inflow, 7in = 1- The 1st and 3rd 
terms are the energy flux in the inflow and outflow region, 
respectively. Note that the 2nd term in left-hand side of the 
equation expresses kinetic and electric field energy of the in¬ 
jected turbulence; The turbulent components in the sheet is 
neglected because we use a sub-Alfvenic turbulence whose 
kinetic energy is small compared with the other terms. Us¬ 
ing the pressure balance: ps = pin-l-Bfn/87r7?, the steady state 
condition: cEy = BinVin = fish's, and the equation of continuity. 
Equation ([T]i, we obtain 


Ps _ (2/Zinfr “i"40)7s . 1-t 2/Zinrr Cinj Ix^z 

Pin /rin(l+cr)-7s _ (2/rintT- 1 - 40)72 S Vsc2_ ’ 


(6) 
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Figures. Reconnection rates fitted by two functions: fiVinj/^A) = 

Cl(Ps/Piii)[Vinj/cA -C2(Vinj/cA)2], g(Vinj/cA) = Ci(ps/pin)Vinj/cA. The func- 
tion / can explain the reconnection rate from the matter dominated case to the 
Poynting dominated case, which indicates the compressible effects definitely 
affect the reconnection rate in turbulence. 


velocity Vin in the left-hand side of Equation Q, and this re¬ 
sults in the decrease of the sheet density comparing with the 
inflow density as indicated in Equation (fTli. In other words, 
this is because the second term in the leftmand side increases 
with Einj but right-hand side only increases with 6 cx . 
Hence, to keep the energy flux conservation, the system re¬ 
duces the inflow velocity Vin in the left-hand side of Equation 
0 , and this results in the decrease of the sheet density com¬ 
paring with the inflow density as indicated in Equation Q. 

Note that Equation (|^ indicates the density in sheets be¬ 
comes negative when a too strong turbulence is injected, 
which does not occur in real situations. This is prohibited 
by including neglected te rms in Equation Q. In particular, 
as is discussed in Section 4.2 energy flux escaping as com¬ 
pressible waves cannot be neglected as the injected turbulent 
Alfven Mach number approaches unity. 


where we neglected a small term proportional to {5 /L)2 ^ 1 
resulting from Note that the denominator, /!in(l -I-ct)- 

7s, is always positive because 7s ^ 7a = vT+ct where 7^ is 
the Lorentz factor of the Alfven velocity in the upstream re¬ 
gion. Since LV99 predicts 8 cx L(vinj/cA) (see Equations (2) 
and (4)), the second term in Equation ([^ becomes propor¬ 
tional to: Cinj/Lvinj cx Vinj; And we finally OTtain the following 
form of relation: psjPin = Qcjl -/3vinj/cA) which qualitatively 
reproduces the linear dependence of the density ratio on the 
injected turbulence strength Vinj in Eigure 4 

This decrease of the sheet density can be explained as fol¬ 
lows. When the turbulence energy injection rate is small, 
Vinj ^ Ca. it increases the sheet width as predicted by LV99. 
However, the increase of the sheet width is proportional to 

I Vinj I fx eU , and the turbulence energy injection rate einj can¬ 
not be absorbed into the sheet width expansion as indicated by 
the second term in Equation Q. In this case, in order to keep 
the energy flux conservation, the system reduces the inflow 

® Our simulation results did not reproduce the exact value of a and /3 
indicated by Equation However, the obtained values of a roughly repro¬ 
duced the predicted value by Equation j^, and our results also reproduced 
the increasingnature of a and P in terms of a parameter; this is indicated 
by Equation with an assumption v, ~ ca if Vj < 0.3c ot herwise v,, = 0.3c 
which is a known result for relativistic Sweet-Parker sheet jTakahashi et al.| 
|201 If Concerning /3, note that it is difficult to estimate the exact value from 
simulation results because of the uncertainty of 4,4, <5. 


4.2. Compressible Turbulence Effects 
Next, LV99 obtained the following relation: S/L (x 
( vi / ca )^ (X Vinj/cA using the incompressible MHD turbulence 
cascade law. In this paper, we treated compressible MHD tur¬ 
bulence, so that it is expected the above relation should be 
modified. More precisely, the LV99’s relation can be rewrit¬ 
ten as: 



(7) 


and substituting, Cinj vf f2lcA, recovers Equation 0. Hence, 
if we find an expression of the energy injection rate Ci^j in¬ 
cluding compressible effects. Equation (j? 1 may give 1 


expression of the sheet width. Recently, Baneijee & Galtier 
(20131 obtained an exact relation of energy cascade rate in 
the non-relativistic isothermal MHD turbulence. In the strong 
background average magnetic field limit, the relation reduces 
to: 

-4-e = V-F+bIS ( 8 ) 


where the divergence V is performed on the correlation length 
which plays a role of the eddy scale length, F is the energy 
flux vector including compressible effects with order of Bq, 
and B is a source or sink term due to the compressible ef¬ 
fects. This indicates that the compressible effects cannot be 



























5 




Figure 6. Top: The ratio of the compressible and incompressible velocity 
components: (v^) and (v^), respectively. Bottom: The ratio of the fast MHD 
wave and Alfven wave power: (V//Va)“- The data were calculated by the 
background without sheet to avoid inhomogeneous back ground and back 
reaction from the sheet. The panel shows the ratio with respect to the strength 
of the injected turbulence: Vinj/cA with various kinds of the magnetization 
parameter. 

neglected in the strong background magnetic field, and the 
energy cascade rate should be redefined as an effective mean 
total energy cascade rate: = e + B^S/A, and this will give 

us the necessary correction term in Equation Q. n Perform¬ 
ing the Taylor expansion of Ceff in Vinj /ca < 1 up to 2nd-order, 
the corrected sheet width can be written as: 


performed using a background obtained by new runs without 
current sheets but using the same setup. This is because the 
current sheets introduce a inhomogeneous background which 
makes it very difficult to perform the above decomposition. 
The top panel of Figure 6 are the ratio between the compress¬ 
ible and incompressible velocity components: (v^) and (v?), 
respectively in terms of the Alfven Mach number of the in¬ 
jected turbulence velocity dispersion, Vinj/ ca, using various 
magnetization parameter cr. As expected, this shows that in 
all cases the compressible component increases with the tur¬ 
bulent strength. We also note that the maximum reconnec¬ 
tion rate in Figure 2 is obtained when (v^)/(vf) ~ 0.4. In¬ 
terestingly, Figure 6 shows the compressible component in¬ 
creases with the magnetization o f the plasma similarly to in 
non-relativistic MHD turbulence ( |Cho & Lazarian|^02[ ). We 
consider this is due to the Bq factor in Equation pjiTwhich 
indicates the compressible effects becomes more important 
as the background magnetic field increases. Unfortunately, 
Equation (|^ is a result in the case of the non-relativistic MHD 
turbulence, and its relativistic extension is our future work. 
The bottom panel of Figure 6 is the ratio of the fast MHD 
wave to the Alfven wave power. The wave decomposition 
was performed in the Fourier space assuming the linear dis¬ 
persion relation si milarly to the non-relativistic case by |Cho| 
|& Lazari^p002| l. It shows a very similar behavior obtained 
by the Helmholtz decomposition. Interestingly, we also found 
the slope of the ratio, (V//Va)^, is proportional to y/a, so that 
it can be written as: (VflVA)^ oc \f(JV\^-^lcA- Equation can 
effectively be derived by considering the wave decomposi¬ 
tion. The turbulent reconnection theory in LV99 considers the 
MHD turbulence results in a wider current sheet because of 
the wandering motion of the magnetic field driven by Alfven 
waves. Hence, Cinj in Equation 0 is equivalent to the Alfven 
wave power, V^, where Va is the Alfven wave component of 
the velocity. In the compressible regime, a part of injected en¬ 
ergy is distributed into the fast wave, and einj in Equation Q 
should be rewritten as: 

^ V?nj-^A V^Vinj/cA, (10) 



where C^ is a coefficient resulting from the expansion. Fig¬ 
ure 5 is the reconnection rates with various kinds of magne¬ 
tization parameters: cr = 0.04,0.5,1,5 which are fitted by 2 
functions; one uses Equation 0 with the density ratio. Equa¬ 
tion §: /(Vinj/CA) = Ci(ps/Pin)[Vinj/cA-C 2 (Vinj/cA)^], and the 
Other only takes into account the density ratio and uses LV99 
sheet width. Equation 0: .g(vi„j/cA) = Ci(ps/Pin)Unj/cA- Note 
that Cl and C 2 describe coefficients independent of the injec¬ 
tion velocity vinj indicated by Equation Q. As can be seen, 
they are well-reproduced only by / whos ^2 are around unity 
in all the cases. 

The above discussion assumes that the turbulence becomes 
compressible. To conform the validity of this assumption, we 
performed the Helmholtz decomposition and the MHD wave 
mode decomposition of the velocity field. The analysis were 

^ The source/sink term 5 from compressible effects includes (1) mode 
exchange between the compressible modes and Alfven mode; (2) a direct 
cascade of eddy size by dilatation (V ■ v > 0) or compression (V ■ ft < 0) 
(more detailed discussion can be found in (Banerjee & Galtier|20131). It is 
our future work to determine which effect works dominantly. 


where Vy is the fast wave component of the velocity, 
reduces to 



Vinj 


^ An] 

2 CA 


This 

( 11 ) 


where we assume Vjnj /ca < 1. This reproduces the depen¬ 
dence of Equation (|Ty on Vinj However, we cannot find 
the dependence of C 2 on ^ indicated by Equation ( [ll] i. We 
consider this may be due to the effect of inhomogeneous nack- 
ground structure and the back reaction from the tearing insta¬ 
bility which are not taken into account to obtain Figure 6. 


5. APPLICATIONS TO HIGH ENERGY ASTROPHYSICAL 
PHENOMENA 

In this section, we discuss applications to high energy 
astrophysical phenomena, that is, the Crab pulsar wind 
nebula, relativistic jets, and gamma-ray bursts. We estimate 
the necessary spatial and temporal scales for explaining those 
phenomena, and compare them with the actual observational 
indications. In those phenomena, it is natural to consider 
the driving source of turbulence depends on phenomena and 


^ Note that slow waves basically propagates along magnetic fields, and 
they are also responsible for the magnetic field wandering. 
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the resulting turbulence strength is different. However, there 
are still a lot of theoretical and observational uncertainties 
in those phenomena, such as the precise spatial distribution 
of the magnetic field strength, the particle number density, 
and the particle composition which is necessary to estimate 
synchrotron radiation flux. And identifying each turbulent 
process is beyond the scope of this paper. Hence, in the 
following, we only use the values of reconnection rate to 
estimate those scales, so that the discussions can in general 
be applied for another dissipation mechanism, such as 
collisionless reconnection or plasmoid-chains. 


-Crab Pulsar Wind Nebula 

In the case of the Crab pulsar wind nebula, the wind region is 
filled with current sheets (striped wind). It is known that the 
magnet ic field cannot be completely dissipated in the wind 
region (|Lyubarsky & Kirk|2001[ |Kirk & Skjaeraasen||2003|l. 
One way to avoid this problem is to assume that the magnetic 


field dissipates just behind the termination shock (Lyubarsky 
2003 1 [Petri & Lyubarsky||2007l Sironi & Spitkovsky [201 Ij l. 
If assuming high-cr upstream now, the down stream of the 


termination shock is still relativistic flow with Lorentz factor 
In the downstream rest frame, the sheet separation is 
around Tv-Jari^c where tlc = c/H is the light cylinder radius , 
and H is the rotation period of the Crab pulsar. The necessary 
timescale to dissipate the magnetic field between the sheets 
by reconnection with reconnection rate Vin is vry^rLc/vin. In 
the termination shock rest frame, which would be equivalent 
to the observer’s frame, additional Lorentz factor is 
multiplied to the dissipation time scale, and it becomes 


Alissip.PWN = TTCTrLc/Vin- 


( 12 ) 


During this timescale, the current sheets propagates 


^dissip “ ^Tiissip.PWN “ 


= 5 X 10“^ [pc] 


/ 27rrLc \ 
^10® [cm]/ 




(13) 


which is sufficiently short compared to the Crab pulsar wind 
nebula scale size (~ 1 [pc]), and the turbulent reconnection 
can be one of the possible dissipation mechanisms to solve 
the (T-problem. The indicated value in Equation ( [T3| l is still 
too small to be resolved by X-ray (e.g. by the Chandra 
telescope). However, it may be possible to be observed by 
a future mission if a is larger than 10"* and the reconnection 
rate Vin/c is smaller than 0.1. 


-Relativistic Jets 

In this case, although there is no proof of the existence of 
current sheets in the observed jets, we assume a dynamo- 
process such as the magneto-rotational instability (MRI) in 


their accretion disk results in current sheets in jets ([Barkov 
& Baushev 201 1]|. Assuming the separation in the fluid 


comoving frame as I, the reconnection timescale can be 
written as: fdissip = l/vm in the fluid comoving frame. If 
we assume the reversing of the magnetic field direction 
occurs every one-Kepler rotation time at a radius r, the 
timescale in the blackhole rest frame can be written as, 
Aot T^rui^r jru)^^^Ic where rm is the Schwarzshild radius. 


Hence, I can be written as, FjetCTrot- Since the dissipation 
timescale in the central blackhole comoving frame, Ljetfdissip. 


should be less than the jet propagation time, /jet/c, we obtain 



^jet 


L jet 




3/2 


19 X If’ (sfcl)"' (s)’'” (»%|) (%i)” tailio). 

3 J X 10-< ((^) (!f)’ (HST), 

(14) 


where the upper and lo wer values are based on th e radio 
and HST data of M87 ([Rieger & Aharonian||2012[). Here 
we assume the reversing of the magnetic held direction 
occurs at the innermost secure radius r = Sr^- This indicates 
the both observation results can be explained by turbulent 
reconnection. 


-Gamma Ray Bursts 

It is suggested that GRB c an be explained by Poynting- 
dominated plasma model, and|Zhang & Yan|P01 l|l suggested 
a model called ICMART model that c an explain GRB includ - 
ing prompt emission spectral curves ( [Zhang & Zhang|2014[ ). 
In the model, the reconnection rate is assumed a relativistic 
value, the minimal value is around 0.1c. As indicated in 
Figure 2, this value of reconnection rate will be obtained 
by turbulent reconnection in Poynting dominated plasma. 
The model also discussed a possible role of reconnection 
outflow assuming Alfvenic velocity, ^ s/l+<J. In our 
calculations, a relativistic outflow was observed only locally 
and intermittently; the averaged outflow velocity was a 
sub-relativistic velocity up to ~ 0.3c even in the high-cr 
regime, which is commonly seen in the relativistic MHD 
reconnection resulted from the tearing instability. If we see a 
smaller scale, the collisionless plasma regime will appear in 
which t he outflow veloc ity is Alfven velocity in the high-cr 

regime ( Liu et al.|20151l. _ 

Recently, [Lazarian & Medvedev] ([2015[) proposed a new 
scenario for the gamma-ray bursts powered by turbulent 
reconnectio n based on kink instab ility of relativistically mag¬ 
netized jets ( [Mizuno et al.||20f^ . In this model, the authors 
considered a turbulence induced by the kink instability, and 
applied it to the turbulent reconnection model. They showed 
that their model can provide a good fit to the dynamics of 
GRBs. Note that this model may also be able to be applied 
to AGN jets, which may result in a different condition from 
Equations ([T4|). 


-General Remarks 

Finally, we give several comments which can be applied to 
any phenomena in general if turbulent reconnection works. 
First, our results, such as the Figure 2, indicates that the 
maximum reconnection rate is obtained when Vinj/cA ~ 0.2. 
Assuming the Alfvenic mode, SB/Bq 5v^/ca, where Bq 
and SB are the background and fluctuation of magnetic field, 
respectively, 5v± is the fluctuation velocity perpendicular to 
Bq, we can expect that the turbulent reconnection is efficient 
if an observed magnetic fluctuation of some phenomena is 
around SB ~ Q.2Bq. In addition, we can expect that such 
a low magnetic fluctuation will allow a high polarization 
degree of synchrotron radiation observed by[Coburn & Boggs] 
(|2003|). Second, it is well-known that magnetic reconnetion 
is also related to particle acceleration. For example, many 
Particle-In-Cell simulations indicates that the energy spec- 
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trum index obtained in relativistic collisionless pair-plasma 


reconnection is around -1.5 ( 

Zenitani & Hoshino 

2001 

Bessho & Bhattacharjee 2012| 

Sironi & Spitkovsky 

2014 

Guo et al.|20151l. On the other hand, |de Gouveia dal 

’ino & 

Lazarian|(|2005|l found that the turbulent reconnection is also 


a location of the first-order Fermi process, and they obtained 
a little softer energy spectrum index, -2.5. Non-thermal 
particle energy index is directly related to the observed 
synchrotron spectrum, and the above energy spectrum index 
may be useful tool to determine the location of turbulent 
reconnection. Applications of turbulent reconnection model 


to some AGNs are provided in (Khiali et al. 2015 Kad 
owaki et'ar]|2015[ [Singh et al.||2015|l. More comprehensive 


discussion of recent observations is given in 

( |2m31 ). 


Lazarian et al. 


6. DISCUSSION AND CONCLUSION 

In this paper, we investigated turbulent reconnection in rel¬ 
ativistic plasmas from the matter dominated to Poynting dom¬ 
inated cases using the relativistic resistive MHD model. The 
results show that the turbulence can enhance magnetic recon¬ 
nection even in relativistic plasmas, and can be a candidate for 
a fast reconnection process. We found the reconnection rate 
in turbulence shows the following 3 characteristic phase de¬ 
pending on the velocity of the injected turbulence: (1) LV99 
region (incompressible turbulence); (2) saturation region giv¬ 
ing maximum rate; (3) reducing due to the compressibility. 
The saturation occurs when the compressible component be¬ 
come dominant, typically around v^/v? > 0.4 at which the 
maximum reconnection rate is about 0.05. This shows that the 
LV99 expressions for incompressible fluid should be modihed 
to account for compressi bility as we have done in E quations 
Q and 0; Interestingly, [Banerjee & Galtierj (|20T3]l showed 
that dilatation of fluid, V • v > 0, reduces the effective energy 
cascade rate Ceff. All of our numerical results showed dilata¬ 
tion, and the reconnection rates are indeed reduced. This indi¬ 
cates that the turbulent reconnection rate may become larger 
than LV99’s prediction if compression of the turbulence oc¬ 
curs, such as t he MHD turbulence driven by collisio ns of mag¬ 
netized blobs pnoue et al.|2011| Deng et al.|2015[ ). 

Finding a fast reconnection process is one of the most im¬ 
portant topics in plasma physics, and a considerable number 
of studies have been conducted on it for a long time. Turbu¬ 
lence is a very general process in high-Reynolds number plas¬ 
mas, so that turbulent reconnection can appear in many kinds 
of phenomena, such as astrophysical phenomena, nuclear fu¬ 
sion, and laser plasma. In particular, our work investigated 
the extension of turbulent reconnection to relativistic plasma 
with compressible turbulence, which allows us to apply this 
process to many high energy astrophysical phenomena, such 
as flares in pulsar wind nebulae, gamma ray bursts, and rela¬ 
tivistic jets. 
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